Anomalous impurity resonance in graphene 



IFIMUP/Departamento 



Yu.G. Pogorelov 

Fisica, Universidade do Porto, R. Campo Alegre, 687, Porto, 4169, Portugal 



o 
o 



A Green function analysis has been developed for quasiparticle spectrum and localized states of a 
2D graphene sheet in presence of different types of substitutional disorder, including vacancies. The 
anomalous character of impurity effects in this system is demonstrated, compared to those in well 
known doped semiconductors, and explained in terms of conical singularities in the band spectrum 
of pure graphene. The criteria for appearance of localized states on clusters of impurity scatterers 
and for qualitative restructuring of band spectrum are established and a phase diagram in variables 
"disorder" vs "electron density" is proposed. 
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PACS numbers: 03.065. Pm; 71.30.+h; 71.55.-i; 81.05. Uv 

There is a growing attention to electronic properties 
of a single carbon layer known as graphene [ij. Its 2D 
honeycomb lattice defines a peculiar band structure 
with two nodal points in the Brillouin zone where conical 
energy surfaces (with zero effective mass) of conduction 
and valence bands touch each other. This gives rise to 
electronic dynamics of relativistic Dirac type jSj , extraor- 
dinary for condensed matter, and generates such unusual 
phenomena as half-integer Hall effect ^ IH 1^ and^pos- 
sibly, the magnetic catalysis of an excitonic gap Q , 
ferromagnetism and superconductivity 0. On the other 
hand, it is of interest to examine the effects that vari- 
ous kinds of impurities can produce on this remarkable 
material, regarding for instance a fundamental role of 
such effects in physics of common semiconductors (with 
finite effective mass) intriguing situation with 

impurity levels near conical singularities was recognized 
in d-wave superconductors where theoretical predic- 
tions are sometimes contradictory |l2j and not fully con- 
firmed by the existing experimental data. To this time, 
the disorder effects in graphene were theoretically stud- 
ied, searching for weak localization in this 2D electronic 
system under weak scattering (Born limit) or for 

strong localization under infinitely strong (unitary limit) 
perturbation ■ This work is aimed on a consequent de- 
scription of restructured electronic spectrum, at arbitrary 
perturbation strength and in a rather broad range of im- 
purity concentration, and on specifics of this restructur- 
ing for Dirac quasiparticles under realistic perturbation, 
compared to usual quasiparticles with parabolic disper- 
sion. 

Let us start from the simplest tight-binding Hamilto- 
nian restricted to nearest neighbor hopping 



H 



n+S 



h.c. 



(1) 



where t is the hopping amplitude, and bn+s are the 
Fermi operators of (spinless) electrons on sites of type 1 
and 2 respectively, and atomic energy on each site is cho- 
sen zero. The vectors S point to three nearest neighbors 
of a site (the vectors for a type 1 site being the inverses 



of those for a type 2 site, see Fig. ^). Passing from site 
operators to plane waves: Ok = N^^/"^ e~*'^ "a„ and 
6k = N~^^'^ e~^^'"bn {N being the number of cells in 
the lattice), and then to the eigen-modes 
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the Hamiltonian, Eq. ^ is diagonalized: 



= ^ Ek ("kak - 0lPv 

k 



with eigen-energies Ek = ^ |/k|- Here the function /k 



„jk<S 



|/k| ^^^^ vanishes near two isolated points in 



the Brillouin zone: K = (0,27r/3a) and -K (Fig. HJ. 
Due to the absence of inversion symmetry for the point 
group C3 , this vanishing is linear in small difference q = 
k - K (or q = k + K): /k = /q ~ ^ (g^; " iqy) defining 
the conical form of isoenergetic surfaces e — ±eq where 



hvpq 



(3) 



with the Fermi velocity vp — \/3ta/2h. The following 
analysis of this system is restricted to the low energy 
physics which is essentially determined by the vicini- 
ties of two nodal points ±K. The four relevant long- 
wave modes, picked up from the eigen-modes, Eq. 12 
near each nodal point, can be combined into the Dirac 
4-spinor ipq whose components (in common notation) 

= /3q+K, 



are: V'q^^ = ct< 



q+K, 



^qi ^ ^' 



q-K, 



<7' - 



= /3q_K- Here, "particles" and "antiparticles" ob- 
viously correspond to electrons and holes and the "Dirac 
spin" indices to the nodal points (while the physical spin 
indices stay suppressed). The respective Dirac form of 
the Hamiltonian, Eq. ^ is 



(4) 
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FIG. 1: 2D lattice structure of a graphene sheet, a) Rhom- 
bic primitive ceU (shadowed) with two non-equivalent po- 
sitions for carbon atoms, 1 (open circles) and 2 (solid cir- 
cles) , and elementary translation vectors ai and a.2 of length 
a = 2AQpA. b) Rhombic Brillouin zone (shadowed) with two 
non-equivalent nodal points, K (open) and — K (solid), and 
vectors of reciprocal lattice bi and b2 of length 47r/(a\/3). 



where the 4x4 matrix 70 = fa (g) ctq, the tensor prod- 
uct of Pauli matrices fi, acting on "particle- antiparticle" 
indices, and ai, on "Dirac spin" indices. We describe 
the dynamics of this system by thefFourier transformed) 
two-time Green functions (GF's) here combined into 
a 4x4 matrix ((V'ql^q'))- For the unperturbed system. 



Eq. ^ its exact form is ((V'qiV'q')) 



with 
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(5) 



The notable distinction of graphene from genuine rela- 
tivistic systems of quantum field theory is the possibility 



to study the effects of localized perturbations on its dy- 
namics, which is our main purpose here. To this end, we 
adopt the Lifshitz model of impurity perturbation [l7|, 
supposing a certain shift V of the on-site energy at ran- 
dom sites of the lattice. This model looks more adequate 
to the case of rare defects in graphene, than the alterna- 
tive choice [1^ of Anderson model with random pertur- 
bations at each lattice cite ^] . We denote pi the defect 
sites of type 1 with concentration ci and p2 those of type 
2 with concentration C2 (not necessarily equal to ci ) , the 
total impurity concentration being C1+C2 = c <C 1. This 
perturbation scatters the modes of Eq. ^ accordingly to 
the Hamiltonian 



H' = 



V_ 

2N 



q,q' I Pi 

P2 ) 



(6) 



with ?7i^2 = (TbT''^i)® (o'o + o'i)- The equations of motion 
for perturbed GF matrix with Hamiltonian H + H' 
can be solved in the T-matrix approximation. For the 
most important, momentum-diagonal GF, this results in 
^1 



matrix is 



Uq 



-S(£) 



where the self-energy 



V 



Its denominator D{e) = 1 + Vg{e) includes the lattice 
sum g{e) = {2e/N) X]k(^^ ~ ^k)^^ which can be approx- 
imated at low energies, |e| <C = \/S7rt/2, as: 



lie) 



In- 



(8) 



A resonance can appear in S(e) at an energy s = Sres 
such that Re D{eres) = 0. But, accordingly to Eq. |S| this 
requires a strong enough perturbation: \V\ > Vcr, where 
the critical value is Vcr ~ 0.947VF. For a numerical check, 
we adopt the common value of i ~ 2.5 eV which leads to 
Vcr ~ 3.6 eV. The defects usually discussed in graphene 
are vacancies , for which the "acceptor" perturbation 
parameter V can be approximated by the 1st ionization 
potential of carbon V Ic ~ 11-3 eV « 31^. It 

results in a low resonance energy: Sres ~ 0.04814^ ~ 0.18 
eV, as shown in Fig. 13 This can be compared with the 
two traditional impurities for carbon compounds: accep- 
tors by boron (1st ionization potential Ib ~ 8.3 eV) and 
donors by nitrogen [In ~ 14.5 eV). The related esti- 
mates for impurity perturbation: Vb ~ Ic ^ Ib ~ 3 
eV and V^v ~ Ic — In ~ —3.2 eV, make it less proba- 
ble that these impurities produce resonances in graphene 
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FIG. 2: Dispersion curves (solid lines) for graphene with ac- 
ceptor impurities at the choice of perturbation parameter 
V — 3W and impurity concentration c — 1%. The dashed 
lines show the unperturbed dispersion laws and impurity res- 
onance level Eres, and the arrows mark the interval of negative 
dispersion. In fact, the solid curves only make sense beyond 
the hatched area of localized states (cf. to Fig. |^) . 



spectrum. Therefore, the foUowing analysis concentrates 
on the effects of vacancies in this material. 

The self-energy matrix, Eq. \7\ can be simplified un- 
der the most natural assumption that defects are equally 
present in two sublattices: ci = C2 = c/2. Then it be- 
comes: 



E (e) ^ cU 



V 



2 [1 + Vgie)] 



(9) 



with U = To ® (ao + ai), that is non- mixing particles 
and antiparticles. This facilitates diagonalization of GF 
matrix, through a spinor rotation: ijjq — Qi^q, with Q = 
To ^ (o'o — io'2) / \/2, leading it to the form: 



\ 







(10) 



with the scalar self-energy S(e) = cy/[l -I- ^^(e)]. We 
notice that two modes, the 1st and 3rd diagonal el- 
ements in Eq. El stay unperturbed. These modes, 
("q+K - aq-K)/V2 and (/3q+K - /3q-K)/y2, are anti- 
symmetric in "Dirac spin" indices, while the perturbation 



U is symmetric in these indices. Thus the impurity scat- 
tering in this case perturbs only two symmetric modes, 
whose dispersion is described by the equations 



£1,2(9) - ReS(ei,2('7)) ± £q 



0. 



(11) 



Using the explicit function 5(e), Eq. |S| and the linear 
law Eq, Eq. 13 we obtain the dispersion curves £1,2(9) as 
shown in Fig. |21 (the picture for "donor" perturbation 
V < Q will be simply inverted, e — > — e, by the particle- 
antiparticle symmetry) . 

If Eqs. El were valid down to g ^ (though in fact 
they are not), they would describe the shift of nodal point 
energy from zero to a finite value £1.2(0) = £0 = Re S(£o)- 
This value grows with the defect concentration as 



£0 



cV, 

Eres (1 - Co/c) 



c < Co, 
c > Co, 



(12) 



where the characteristic concentration co ~ [eres/WY' 
defines, as will be seen below, the threshold for quali- 
tative restructuring of the spectrum. For c > co, there 
appears a certain interval of negative dispersion in the va- 
lence band (a Z-like feature in Fig. j^J. It resembles the 
known situation near resonances in parabolic bands |20l | , 
however in this case negative dispersion appears rather 
far from the resonance Sres (which lies within the con- 
ductance band). And the restructured conductance and 
valence bands are both shifted towards this resonance, 
not repelled from it (as commonly for doped semicon- 
ductors). 

To clarify the physical origin of this anomalous behav- 
ior, it is instructive to compare it with the well stud- 
ied case of shallow donor levels below a 2D parabolic 
conductance band Ek = h^kF' /2m [l^. The respective 
scalar self-energy formally coincides with the above ex- 
pression for E(£) but including the lattice sum gie) = 
iV~^^^(£ — Ek)^^ w ln(i?„,./£) and a weak attractive 
perturbation < — <C E„i ~ ^-kY? j ma} . Its expan- 
sion near the localized donor level £;oc ~ — i?me~'^'"/'^' 
reads: E(£) w cV\eioc\l{^ — e/oc), so that the perturbed 
band energies e = E^ + E(£) are repelled from this level. 
In fact, this repulsion is due to the simultaneous action 
of attractive perturbation V in the numerator and de- 
nominator of S(£) (so it remains also true for repulsive 
perturbation, V < 0, hy shallow acceptors). In contrary, 
such expansion for the actual case of linear band dis- 
persion reads: E(£) « —cV\eres\/{£ — Eres), like if the 
repulsive perturbation > in the denominator of E(£) 
turns to be attractive in its numerator (or vice versa for 
donors), resulting in the overall attraction of the bands 
to the impurity resonance. 

However, the above referred dispersion curves are only 
reliable when the respective band states are well defined, 
which can be checked by the loffe-Regel-Mott (IRM) cri- 
terion [2^ that the quasiparticle lifetime is long enough 
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FIG. 3: Check for the IRM criterion at different concentra- 
tions of defects, a) Below the critical concentration, c = 
0.1% < Co, all the states are band-like except a very narrow 
vicinity of the nodal point energy eo (shown under a great 
magnification in the inset), where the states are localized, b) 
At approaching the critical concentration, c — 0.2% « co, 
the area of localized states rapidly expands on the whole in- 
terval between eo and Sres- c) At overcritical concentration, 
c = 1% S> Co, the area of localized states of width F extends 
beyond the resonance level Sres- 

compared to the oscillation period. This can be suit- 
ably presented in the energy scale by the inequality 
q(e) > K{e) where the "real wave number" q{e) is the 
inverse function to e{q) and the "imaginary wave num- 
ber" = Imi;(e)[l -Re !]'(£)] /ftwp. As seen from Fig. 
13 the IRM criterion ceases to hold just for energies close 
to Eq, where the quasiparticles are not properly described 
by the quasimomentum and band index, but rather lo- 
calized near impurity centers. Therefore the band states 
in the perturbed system never include the nodal points. 

A closer insight on the localization process can be ob- 
tained, beginning from the case of a single impurity on 
site p (say, of type 1) and constructing a quasiparticle 



state with energy e: lipe) = J2n (tf ''fln + V'i^^^n) |0), 

(1 2) 

whose amplitudes on type 1 and 2 sites ipn ' are found 
from the Schroedinger equation, {H + H' — e)\ipe) — 0. 
The solutions having central symmetry with respect to 
the impurity site p: 

T/ . . ;.„iq (p— n) 

q ^ 

= <V((ap|al,)), 

q ^ 

= <V((ap|6t,)). (13) 

are proportional to the locator Green functions. Their 
long distance asymptotics (at |n — p| a) is: 

with the characteristic length r{e) = {Wa)/{^/^Te) ^ a. 
In particular, the value tq = r{eres) defines the length 
scale of local perturbation of quasiparticle spectrum near 
a defect. So, a qualitative restructuring of this spectrum 
should happen if local perturbations begin to overlap, at 
characteristic concentration of defects: cq ~ /{-rtq) = 
{Erea/WY. For the considered vacancy model, this value 
is ^ 0.2%. In fact, when the concentration c reaches this 
level, important changes occur in the spectrum charac- 
teristics calculated from the GF, Eq. ^1 

For instance, the quasiparticle density of states 
(DOS) is p{e) = (7r7V)-ilmTrEq((V^q|^t)). At c « 
Co, it is close to the unperturbed DOS po(£) = 
(7riV)-ilmTr^qG° = |e|/VF^ but at c ~ cq a hump 
appears in p{e) near the resonance Sres and then pro- 
gressively increases at c > co (Fig. 0J. 

Next, the Fermi energy ep is obtained in function of c 
from the condition ^^^p[e)de = (1 — c)/2, and (in the 
given model) it displays a rapid initial growth at c ^ cq 
(inset in Fig. 01, entering the conductance band and 
so realizing an anomalous n- conductance at nominal p- 
doping. But, before c approaches cq, the level ep makes 
a sharp downturn and crosses the monotonically growing 
nodal energy eo, thus restoring a usual type of conduc- 
tance. 

At last, the width F of the energy interval near eo 
(hatched in the inset in Fig. , where the IRM criterion 
ceases to hold, is as narrow as ~ eoc/(co ln(l/co)) ^ eo 
at c ^ Co but suddenly expands at c ~ co to the whole 
range from eo to e^es, and then grows as ~ -yc/ln(l/c) 
at c ^ Co. This energy interval is filled by localized 
states, however their localization is realized as a rule on 
certain clusters of defects, rather than on single defects. 
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FIG. 4: Spectrum restructuring with growing concentration 
of vacancies (indicated by numbers at solid lines for DOS 
p(e)). Dashed lines show the unperturbed DOS po(e) and the 
resonance energy Sres- Inset: concentration dependencies of 
the Fermi energy ep and nodal point energy £o (the hatched 
area of width V represents its concentration broadening). 

On each side of this interval, the localized states are sepa- 
rated from the band-like states by a Mott's mobility edge, 
defining a metal-insulator transition when the Fermi en- 
ergy crosses this edge. 

It is of interest to compare this type of spectrum 
restructuring and the quadratic law co ^ (Sres/W^)^ 
with the known examples for non-relativistic spectra. 
Thus, the low-frequency acoustic resonance lOres ~ 
\/ M/M' <^ uju by impurities with mass M' ^ M 
in a crystal with atomic mass M and Debye frequency 
'^D 13 gives rise to splitting of phononic spectrum near 
ujres with opening of a quasi-gap (seen as a dip in DOS) of 
width coj'^/ujres, at surpassing the characteristic irnpu- 
rity concentration cq ~ {oOres/^'DY' (the cubic law) |25| . 
This dip corresponds to repulsion of band levels from 
the impurity resonance (in contrast to the hump in DOS 
and attraction of band levels to the resonance expected in 
graphene). Another example is the donor level ejoc near a 
parabolic conductance band (22j, which rapidly expands 
and merges with this band when the donor concentration 
exceeds Qoc ^ sioc/Em (linear law), in the 2D case, or 
cioc ~ {^loc/ Em)"^^"^ , in the 3D case. The 3/2 law also de- 
fines the characteristic concentration cq ^ {cjres/ J)^^^ of 
weakly coupled (J' <ti J) impurity spins in a Heisenberg 
ferromagnet when the magnon spectrum splits near the 
resonance frequency tOres ^ J' J , and the cubic law 
cn ~ [ujres/JY defines such effect in an aniferromagnet 
|20| . Hence the case of defects in graphene differs from 
all those, even at seemingly identical linear dispersion (as 
for phonons and antiferromagnons). 

Observable effects related to the restructured elec- 



FIG. 5: Phase diagram of electronic states in graphene under 
disorder and variable electron density (with zero level corre- 
sponding to exact half filling. 

tronic spectrum of graphene with defects need a detailed 
consideration which shall be done elsewhere. Here we 
only outline a general framework for this study, through 
the phase diagram in concentration of (vacancy) defects 
c vs electron density n (the latter can be independently 
tuned, e.g., by the external bias voltage or by adding 
non-resonant donor and acceptor impurities). Accord- 
ingly to the concentration dependence of characteristic 
energy parameters in the inset in Fig. 01 this diagram 
(Fig. 0) presents three different phase regions: the "in- 
versely doped" metallic phase Mi with e-p > Eq, the 
insulating phase / with e-p inside the domain F of lo- 
calized states, and the "normally doped" metallic phase 
M2 with sp < Eq. Within the latter phase, a specific 
region of inverse, negative dispersion N can be distin- 
guished, where the Fermi carriers should behave like holes 
by their electric charge but like electrons by their veloc- 
ity, the latter turning extremely high near the inversion 
line between M2 and N. This peculiar dynamics can 
produce even more extraordinary anomalies under ap- 
plied magnetic field (which effect is already anomalous 
in the pure graphene). The latter factor and also inter- 
action between electrons in presence of defects need a 
more involved treatment. All the mentioned features can 
be potentially used for applications in a specific "rela- 
tivistic" electronics, based on properly doped graphene, 
and undoubtedly deserve further attention. 

In conclusion, the effects of local perturbation by var- 
ious types of impurities, including vacancies, in single 
layer graphene were studied through Green function tech- 
niques. Possibility of low energy resonance near nodal 
points in the relativistic electronic spectrum was indi- 
cated for the case of vacanicies and the conditions for 
qualitative restructuring of the quasiparticle spectrum 
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were established. A significant distinction of this restruc- 
turing for relativistic electrons, compared to the known 
disordered systems, is demonstrated and the phase dia- 
gram in variables "disorder vs electron density" is pro- 
posed, indicating possible practical applications of such 
system. 
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